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Abstract 


- 


This  study-evaluated  four  possible  technique^for 

r 

combining  single  axis  pilot  ratings  to  predict  a  pilot 
rating  for  the  axes  in  combination.  Two  of  the  methods 
proved  promising,  giving  good  results  for  the  experi¬ 
mental  data  taken  during  the  study.  The  first,  developed 
by  I.  L.  Ashkenas,  is  based  on  a  non-linear  relationship 
between  the  pilot  rating  and  the  pilot  lead  time- const ant. 
The  second,  based  on  the  incremental  change  of  pilot 
rating  as  the  vehicle  handling  qualities  are  changed, 
gives  comparable  results,  but  is  restricted  as  to  the 
data  to  which  it  can  be  applied.  An  attempt  to  improve 
the  results  of  the  two  methods  by  combining  certain 
features  of  each  proved  to  be  unfruitful.  The  fourth 
method,  based  on  the  summation  of  the  individual  axes 
ratings,  was  shown  to  be  useless. 
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AN  EVALUATION  OP  FOUR  METHODS  FOR  CONVERTING 
SINGLE  AXIS  PILOT  RATINGS  TO  MULTI-AXIS 
PILOT  RATINGS  USING  FIXED-BASE  SIMULATION  DATA 

I .  Introduction 

Because  the  term  ‘'handling  qualities”  has  a  vague 
rather  than  concrete  meaning,  assessment  of  these  qualities 
is  subjective  rather  than  objective  or  quantitative.  How¬ 
ever,  experience  has  shown  that  such  assessments  are  depend¬ 
ent  on  definable  vehicle  stability  and  control  character¬ 
istics.  To  the  pilot,  the  three  characteristics  that  are 
of  primary  importance  are  the  closed  loop  response,  the 
open  loop  characteristics,  and  the  maneuvering  capabilities 
of  the  system. 

Of  particular  interest  is  the  closed  loop  response 
because  the  pilot  is  an  active  participant  in  the  control 
system,  acting  as  a  vocal,  adaptive  controller.  Since  he 
is  aware  of  his  role  and  performance  in  the  system,  his 
opinion  (quantized  as  a  "pilot  rating”)  of  its  handling 
qualities  is  a  dominant  consideration  in  assessing  the 
overall  pilot-vehicle  performance. 

A  pilot-vehicle  system  theory  of  handling  qualities 
has  been  developed  and  extensive  analysis  (Ref  3,  5)  has 
been  made  for  single  axis  situations,  Ihe  method  has 
progressed  to  the  point  of  being  able  to  predict  pilot 
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ratings  for  single  axis  tasks.  During  the  development 
and  verification  of  the  theory,  a  considerable  amount 
of  single  axis  experimental  data  has  been  accumulated. 

However,  actual  flight  situations  require  simultane¬ 
ous  control  of  more  than  one  axis.  In  extending  the 
single  axis  theory  to  multiple  axis  tasks,  it  would  be 
desirable  to  utilize  the  information  gathered  in  the 
single  axis  experiments  rather  than  to  simulate  a  large 
number  of  multi- axis  configurations.  The  question  arises, 
however,  as  to  how  to  combine  the  single  axis  data  to 
determine  a  rating  for  a  multi-axis  task. 

The  object  of  this  thesis  was  to  evaluate  various 
techniques  that  might  be  used.  The  approach  was  to  gather 
single  and  related  multi-axis  data  and  then  compare  the 
multi-axis  ratings  predicted  by  each  method  with  the 
actual  multi-axis  ratings. 

The  pilot's  rating  of  the  handling  qualities  of  a 
vehicle  is  influenced  by  the  dynamic  characteristics  of 
the  vehicle.  In  conventional  aircraft,  helicopters,  and 
spacecraft,  both  electronic  and  mechanical  stability 
augmentation  systems  are  used  to  improve  the  dynamic 
characteristics  when  the  aerodynamic  response  character¬ 
istics  are  poor.  Thus,  it  is  possible  to  classify 
changes  in  the  handling  qualities  of  a  "simulated" 
vehicle  as  relative  changes  in  the  degree  or  amount  of 
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augmentation.  In  this  study  a  low  level  of  augmentation 
implies  the  poorest  handling  qualities  being  considered, 
while  a  high  level  implies  the  best  handling  qualities 
evaluated. 

Scope  of  Study 

This  thesis  evaluated  four  possible  techniques  for 
converting  pilot  rating  (PR)  data  for  separate  single 
axes  into  pilot  ratings  for  the  axes  in  combination. 

The  study  was  limited  to  an  uncoupled,  three  axis,  fixed 
base  simulation.  The  vehicle  transfer  characteristics 
were  derived  f  rom  the  linearized  equations  of  motion  of 
a  tandem  rotor  helicopter  in  "up  and  away",  straight  and 
level  flight.  The  only  criterion  for  the  "goodness"  of 
a  particular  configuration  was  the  pilot  rating. 

Plan  of  Development 

The  pilot  model  transfer  function  and  adaptation  tech¬ 
niques  for  a  compensatory- tracking  task  are  reviewed 
briefly.  A  modified  form  of  the  Cooper  scale  is  presented 
and  comments  made  as  to  why  it  was  used.  The  four  methods 
to  be  evaluated  are  explained  in  detail  with  illustrative 
examples.  The  experimental  set-up,  run  procedures  and 
sequence,  and  operator  training  are  discussed.  The 
methods  used  to  analyse  the  data  are  explained.  Compar¬ 
isons  of  the  four  methods  and  recommendations  for  further 
study  are  made. 
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II.  Pilot  Representation  and  Cooper  Scale 

The  generally  accepted  servo  representation  of  a 
closed  loop  control  system  with  the  pilot  performing  a 
compensatory- tracking  task  is  shown  in  Figure  1. 


In  a  compensatory- tracking  task  the  pilot  observes  only 
the  displayed  error  Di  ,  which  is  the  difference  between 
the  input  and  the  airframe  output  A@.  Hie  pilot's 
output  is  a  force  F  which  produces  a  control  surface 
deflection  S  through  the  control  device.  The  pilot 
transfer  function  Gp  is  usually  defined  as 


GP(S) 


F(S)  -  Ko  ers  (TiS  +  0 
D(S)  ~  (Tt  5  +  DO.S  +  0 
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where 


Kp  =  Pilot  gain 
T  =  Pilot  reaction  time  delay 
Tl  =  Pilot  lead  time-constant 
Tj  »  Pilot  lag  time- const ant 
TN  s  First-order  time-constant  approximat¬ 
ing  pilot's  neuro-muscular  lags 
When  closing  the  loop,  the  pilot  adopts  the  form  of 
equalizing  characteristics  that  insures  stable  control. 

Once  good  low  frequency,  closed-loop  response  has  been 
achieved  with  the  lag  term,  he  generates  whatever  lead 
is  necessary  to  retain  high  frequency  system  stability 
(Ref  3i7).  The  pilot  lead  time-constant  associated  with 
this  generated  lead  is  an  important  parameter  used  in 
the  solutions  for  predicted  pilot  rating  in  two  of  the 
methods  to  be  evaluated. 

Cooper  Scale 

In  order  to  have  some  common  standard  for  aircraft 
handling  qualities,  G.  E.  Cooper  proposed  a  pilot  opinion 
rating  scale  in  1957  (Ref  4:48).  This  scale  has  become 
generally  accepted  for  test  flight  and  simulator  evaluation 
of  flying  vehicles.  A  modified  form  of  the  Cooper  scale 
(Table  I),  used  by  Systems  Technology,  Inc.  in  their 
handling  qualities  study  of  single-loop  roll  tracking 
tasks,  was  used  in  this  experiment.  The  adjective  terms 
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of  this  scale  were  nore  useful  to  a  pilot  with  a  relatively 
low  level  of  flying  experience  than  those  of  the  original 
scale. 
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III .  Explanation  of  Methods 


Method  £1 

The  simplest  method  of  combining  two  or  more  single 
axis  pilot  ratings  into  a  predicted  pilot  rating  for  the 
combined  axes  is  to  sum  all  the  single  axis  ratings. 

The  general  mathematical  expression  for  this  method  is 

cp  r)„  =  y  (p  r); 

where 

(PR)U  a  Pilot  rating  for  the  multiple  axis  configuration 
(PR)i  =  Pilot  rating  for  the  ith  single  axis  configur¬ 
ation 

N  *  Number  of  single  axes  making  up  the  multiple 
axis  configuration. 

It  should  be  apparent  that  this  method  is  not  too 
promising.  In  the  example  in  Table  II,  three  single  axis 
ratings  in  the  ,,good,,  flight  condition  region  sum  to  a 
predicted  rating  in  the  "unacceptable”  region.  It  is 
not  likely  that  the  pilot  would  consider  the  handling 
qualities  of  the  vehicle  had  deteriorated  that  much. 

However,  the  method  is  evaluated  in  this  thesis  for  com¬ 
parison  with  the  other  methods. 
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Table  II 

Example  of  Method  #1 


Axis 

(i) 

Pilot  Rating 

1 

(PR)!  =  2 

2 

(PR)2  =  2 

3 

(PR)3  =  2.5 

1,  2, 

and  3 

(PR)m  =  6.5 

Method  #2 

A  second  possible  technique  for  predicting  multi¬ 
axis  ratings  is  to  take  into  account  only  the  incremental 
effects  on  the  pilot’s  rating  of  a  change  in  the  augment¬ 
ation  in  any  axis  or  axes.  The  sum  of  these  increments 
car.  be  arded  to  some  measured  or  estimated  "best”  multi¬ 
axis  configuration  to  find  the  predicted  pilot  rating. 

The  general  equation  for  this  method  is 


where 


(pR)„=  +  Y. (APR)i 


(PR)  Best  = 


(APR)i  = 


Pilot  rating  for  the  "best"  multiple  axis 
configuration 

Incremental  pilot  rating  between  the  best 
i**1  single  axis  configuration  and  the  i**1 
single  axis  configuration  making  up  the 
multiple  axis  configuration 
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In  the  example  of  Method  #2  shown  in  Table  III 

H  =  Highest  level  of  augmentation  in  any  axis 
M  =  Medium  level  of  augmentation  in  any  axis 
N  =  3 


Table  III 

Example  of  Method  #2 


|  Augmentation 

m 

(PR^Best 

Predicted  PR 

■ 

2 

3 

E 

2 

3 

_ 

BHSsi 

Com 

*• 

H 

H 

H 

B 

B 

B 

3 

B 

H 

H 

M 

B 

B 

a 

3 

4 

C 

H 

M 

M 

B 

2 

1 

3 

6 

In  configuration  A  the  (PR)5est  is  established  from  the 
measured  pilot  rating,  in  this  example  3,  for  the  three 
axes  in  combination  with  the  augmentation  of  each  axis  set 
at  the  "H"  level.  In  configuration  B  the  measured  change 
of  rating  when  going  from  the  "H"  to  "M”  level  in  axis  3 
is  1  rating  point.  This  increment  -  (APR)3-  is  added  to 
the  (PlOfiest  rating  giving  a  predicted  rating  of  4.  In 
configuration  C  the  augmentation  of  axis  2  has  also  been 
changed  from  "IP*  to  "M",  but  with  a  measured  change  of 
rating  of  2  in  this  case.  The  predicted  rating  is  6  or 
the  sum  of  the  (PR)Best  rating  of  3  plus  the  (APIO3  of 
1  and  the  (APR)2  2. 
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Method  #3 

A  tlii  rtl  possibility  is  that  a  predicted  rating  can 
be  found  through  a  non-linear  variation  of  pilot  rating 
with  another  measurable  parameter.  Such  a  method  has 
been  proposed  by  I .  L.  Ashkenas  (Ref  11:10).  Prom  his 
studies  of  longitudinal  handling  qualities  (Ref  9:26-31) 
he  has  found  that  there  seems  to  be  a  non-linear  relationship 
where  the  change  of  pilot  rating  is  a  function  of  the 
pilot  generated  lead  time  constant  TL.  The  graphical 
representation  of  the  function  is  a  band  curve  approximately 
two  rating  points  wide  (Figure  3),  taken  from  reference  11. 
In  this  thesis  the  center-line  of  the  curve  was  used  for 
calculations. 

The  general  equation  for  this  method  is 

cv,...*  {npRM-<TL )....} 

whe  re 

PR  =  f  [  TL  ]  and  TL  =  f"1  [  PR] 

define  the  non-linear  relationship  between  PR  and  Tj^ 
shown  in  Figure  3  and 

('It^Base  =  [Best  single  axis  ?r] 

A  graphical  solution  for  method  #3  is  shown  in 
Figures  2  and  3.  Because  axis(2)has  a  higher  pilot 
rating  than  the  other  axis  at  the  highest  augmentation 
level  (H),  a  PR  base-line,  stuwn  by  a  dashed  line,  is 
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(Adapted  from  Ref  1 1  ) 


Pilot 

Rating 


Figure  2 

Variation  of  Pilot  Rating  with  Damping  Level 


Pilot 

Rating 


Figure  3 

"Ashkenas  Curve" 

Variation  of  Pilot  Rating  with  Pilot  Lead  Time-constant 
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established  at  this  rating  and  is  transferred  to  Figure 
3  at  the  same  value  of  pilot  rating.  The  intersection 
of  the  PR  base-line  and  the  center-line  of  the  conversion 
curve  determines  the  base-line  for  increments  of  pilot 
lead  time-constant  (dotted  line).  The  predicted  pilot 
rating  is  then  determined  as  follows: 

1.  Find,  at  the  level  of  augmentation  of  each  axis 
in  Figure  2,  the  PR  difference  (APR)  between  the 
individual  axes  (l)  and  (2)  and  the  PR  base¬ 
line.  In  this  illustration  axis  (T)  is  measured 
at  the  lowest  level  (L)  and  axis  (2)  is  measured 
at  the  medium  level  (M). 

2.  Lay  off  the  values  of  APR  from  the  PR  base-line 
to  the  center-line  of  the  conversion  curve  of 
Figure  3  and  determine  the  corresponding  pilot 
lead  time  increments  (ATLl  and  ATj^)  from  the 
pilot  lead  time  base-line. 

3.  Add  the  increments  of  pilot  lead  time  and  lay 
off  this  value  (ATj^)  from  the  pilot  lead 
time  base-line  to  the  conversion  curve.  The 
pilot  rating  at  the  point  of  intersection  is 
the  predicted  pilot  rating  for  the  combined 
axes. 
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A  fourth  method  is  to  combine  the  techniques  of 
methods  #2  and  # 3.  The  incremental  changes  of  pilot 
lead  time-constant  are  determined  from  the  changes  of 
pilot  rating  in  each  axis  as  a  result  of  the  augment¬ 
ation  changes.  The  predicted  pilot  rating  is  then 
determined  by  adding  the  appropriate  increments  to  a 
pilot  lead  time-constant  base-line  that  is  determined 
from  a  measured  or  estimated  "best"  combined  axis  pilot 
rating.  The  general  equation  for  method  #4  is 

=f  [0LU+  I 

where 

^TL^Best  =  ^  [  (pR)Best  ] 

Method  #4  is  graphically  illustrated  in  Pigures  4 
and  5.  In  this  example,  the  "best"  pilot  rating  for 
the  combined  axes  is  3  (a  measured  value  in  this  thesis) 
and  the  PR  differences  are  both  measured  at  the  medium 
augmentation  level.  To  determine  predicted  pilot  rating: 

1.  Pollow  the  procedures  for  method  #3  through  step 

2. 

2.  Establish  a  combined  A  Tj,  base-line  based  on  the 
"best"  combined  axis  rating  by  drawing  a  vertical 
line  on  Figure  5  (alternate  dots  and  dashes)  from 
the  intersection  of  the  "best"  combined  axis 
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pilot  rating  line  (alternate  dots  and  double 
dashes)  with  the  center-line  of  the  conversion 
curve. 

3.  Add  the  increments  of  pilot  lead  time  and  lay 
off  this  value  (ATL3)  from  the  combined  A 
base-line  to  the  conversion  curve.  The  pilot 
rating  at  the  point  of  intersection  determines 
the  predicted  pilot  rating  of  the  combined  axes. 
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IV.  Expe rimental  P rocedure 
Equipment  Set-up 

A  servo  model  of  the  closed  loop,  compensatory- 
tracking  control  system  used  in  this  experiment  is 
shown  in  block  diagram  form  in  Figure  6.  The  three 
uncoupled  loops  represent  pitch,  ©  ,  bank,  4>  ,  and  side¬ 
slip,  ^  .  Sideslip  was  used  instead  of  yaw  in  an  attempt 
to  relate  the  dynamics  more  closely  to  an  output  that 
would  be  a  function  of  just  rudder  control.  The  block 
labeled  Gp  represents  the  pilot  transfer  function,  while 
Ge,  G^,  and  G^  represent  the  transfer  characteristics  of 
the  controlled  element  in  each  loop. 


Figure  6 

Servo  Model  with  Three  Uncoupled  Loops 
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The  inputs  into  the  system  (  ©t  ,  4>L  ,  )  were 

generated  by  random  signals  which  were  compared  with 
the  outputs  of  the  respective  feedback  loops.  An 
error  signal  (  ©£  ,  )  was  fed  as  a  voltage  to 

the  display.  The  display  for  pitch  and  bank  error 
appeared  as  a  2  inch  horizontal  bar  on  an  oscilloscope. 
Pitch  error  was  indicated  as  a  vertical  displacement  of 
the  bar,  up  or  flown,  from  the  center  of  the  oscilloscope. 
Bank  error  was  displayed  as  a  tilting  of  the  bar  about 
the  vertical  center-line  of  the  scope.  Horizontal  and 
vertical  center-lines  were  scribed  on  an  overlay  on  the 
face  of  the  oscilloscope  with  intervals  marked  every  .1 
inch.  Sideslip  error  was  displayed  as  a  needle  deflec¬ 
tion  on  a  dial  directly  below  the  scope.  The  face  of 
the  dial  was  marked  at  intervals  to  the  left  and  right 
of  the  vertical  zero  position;  however,  no  significance 
was  given  to  the  numbers  on  the  dial  and  each  mark  was 
considered  an  arbitrary  unit.  Because  the  pilot  found 
it  easier  to  fly,  an  "outside-in"  display  was  provided 
on  the  scope.  Thus,  the  bar  was  flown  back  to  the 
horizontal  center-line  of  the  scope.  He  also  found 
it  easier  to  associate  the  deflection  of  the  sideslip 
needle  to  one  side  with  the  rudder  pedal  on  that  side 
so  an  "inside-out"  presentation  was  used  in  sideslip. 

The  angular  deflection  of  the  sidestick  and  fore- 
aft  displacement  of  the  rudder  pedals  by  the  pilot 
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produced  an  output  voltage  in  each  channel  directly 
proportional  to  the  amount  of  the  movement.  These 
voltages  were  fed  to  the  appropriate  controlled  element. 

The  controlled  elements  were  simulated  on  an 
analog  computer.  One  parameter  of  each  controlled 
element  was  changed  to  produce  three  different  response 
characteristics  in  each  channel.  The  three  were  con¬ 
sidered  as  high  (H),  medium  (M),  and  low  (L)  levels  of 
augmentation  in  each  channel.  The  transfer  functions 
used  for  the  controlled  elements  were 


G  (S)  -  -  4{'S4-.o4)rs  +.9) _ 

"  "  S.s(S)  _(S  +  l/TPe)CS+5)[sil+2(.7)(.25)5+.25a] 


with  l/lj;  =  0  for  H,  =  -.5  for  M,  and  \ /j  -  -*8 


&,(S)  = 


<ft(S)  _  .5(5  +.1) _ 

&.,(&)  (5  +  1.5)  [S!  +  <?y.5)S  +  .5'] 


with  -  +.84  for  H,  f  =  -.3  for  M,  and  ^  =  -.84  for  L 


r  ro  -  _  IO(S  +.0 _ 

'  ;  ~  (5  +  3)[Sa  +  <?S,CS}S  +  .S’] 

with  f  =  +.5  for  H,  f  =  -.2  for  M,  and  <  =  -.5  for  L 


The  input  signals  were  summations  of  non-harmonic, 
low-frequency  sine  waves.  Four  waves  were  summed  tor 
pitch  and  bank  and  three  waves  summed  for  sideslip 
(Figure  A-2).  The  amplitudes  and  freouencics  used 
produced  a  signal  sufficiently  random  for  the  experiment. 
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cillograph  Recording  of  Random  Input  Signals 


G£/ii£/62-4 

The  maximum  amplitudes  of  the  signals  produced  approx¬ 
imately  +  .75  inches  of  displayed  pitch  travel,  +45  degrees 
of  displayed  bank  angle,  and  +1.8  units  of  sideslip 
needle  deflection.  Figure  7  is  a  sample  simultaneous 
trace  of  the  input  signals. 

Run  Procedures 

The  pilot  was  instructed  to  try  and  maintain  the  bar 
level  at  the  intersection  of  the  scope  center-lines  and 
to  maintain  the  sideslip  needle  pointer  at  zero.  After 
a  run  had  been  completed  with  a  particular  configuration 
of  the  augmentation  levels,  the  pilot  was  asked  to  rate 
it  on  the  basis  of  his  interpretation  of  the  modified 
Cooper  scale.  If  he  lost  control  of  the  simulator  dur¬ 
ing  the  run,  the  computer  was  stopped  prior  to  overload. 
The  pilot  was  then  given  the  oction  of  repeating  the 
same  configuration  again  until  he  was  satisfied  he  could 
rate  it.  It  was  felt  that  this  would  give  better  differ¬ 
entiation  of  the  worst  cases.  It  was  also  assumed  that 
the  pilot  would  take  into  account  in  his  rating  any 
improvement  in  his  ability  to  fly  the  configuration  due 
to  fore-knowledge  from  the  previous  run.  He  was  not 
advised  during  the  course  of  the  experiment  doout  control 
techniques,  the  nature  of  the  augmentation  and  how  it 
was  derived,  or  the  consistency  of  his  ratings.  A 
brief  warm-up  period  was  given  prior  to  each  session 
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and  frequent  rest  periods  were  scheduled  throughout 
the  runs. 

Operator  Training 

Time  limitations  precluded  the  training  of  more 
than  one  pilot  for  the  experiment.  The  pilot,  an 
aeronautical  engineer,  had  approximately  180  flying 
hours  in  light  aircraft  and  about  10  simulation  hours 
on  previous  studies.  In  the  three  days  prior  to  the 
data  runs  he  flew  150  training  runs  in  the  single  and 
three  axis  configurations.  The  sequence  of  these  runs 
was  random  and  a  sufficient  number  of  repeated  runs 
were  included  to  determine  the  consistency  of  his 
ratings.  The  values  used  for  the  augmentation  levels 
were  not  those  used  for  the  data  runs. 

An  analysis  of  the  first  training  runs  showed  a 
definite  tendency  to  give  a  much  better  rating  to  a 
repeated  configuration  as  the  pilot  became  more  pro¬ 
ficient  in  flying  the  simulator.  However,  at  the  end 
of  the  third  day  he  was  consistently  rating  repeat 
runs  to  within  an  increment  of  +1.  This  was  considered 
a  reasonable  spread  and  was  compensated  for  to  some 
extent  in  the  data  runs  by  averaging  the  ratings  for 
any  particular  configuration. 
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Data  Kun  Sequence 

Data  for  the  three  axis  analysis  were  taken  during 
the  two  day  period  immediately  following  the  training 
runs.  Twenty-one  single  axis  run--  of  120  seconds  dur¬ 
ation,  followed  by  fifty-three  axis  runs  of  120  seconds 
duration,  were  flown  the  first  day.  The  next  day 
thirty-six  three  axis  runs  of  120  seconds  duration 
were  followed  by  twenty-four  single  axis  runs  of  60 
seconds  duration.  All  single  axis  runs  were  completed, 
but  twenty-eight  of  the  three  axis  configurations  were 
not  completed  on  the  first  attempt.  The  pilot  exercised 
the  re-run  option  thirty-six  times,  repeating  some  con¬ 
figurations  more  than  once. 

Because  the  three  axis  cases  showed  a  preponderance 
of  data  at  the  higher  end  of  the  rating  scale,  a  group 
of  two  axis  runs  with  single  axis  checks  were  flown 
after  a  three  day  delay.  Eighteen  single  axis  runs 
(60  seconds)  and  twenty- nine  two  axis  runs  (120  seconds) 
were  flown  in  one  day.  Unfortunately,  prior  to  these 
runs,  it  was  necessary  to  change  a  turret  resistor  in 
the  pitch  circuit  of  the  computer.  Although  the  same 
numerical  setting  was  used  on  the  replacement,  there 
was  a  definite  overall  improvement  in  the  pitch  char¬ 
acteristics.  The  single  axis  runs  of  this  session  show 
that  the  relative  rating  position  of  each  augmentation 
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level  remained  the  same,  but  there  was  a  large  change 
in  the  increment  of  pilot  rating  between  each  level. 
Since  the  bank  ratings  remained  essentially  unchanged, 
it  was  assumed  that  the  change  in  ratings  was  due  to 
a  change  in  the  circuit  and  not  a  variation  in  the 
pilot.  The  single  axis  ratings  taken  prior  to  the  com¬ 
ponent  change  indicate  that  the  resistor  was  a  fixed 
value  during  the  three  axis  data  runs  and  was  not  vary¬ 
ing  throughout  the  experiment. 

Since  the  procedures  evaluated  in  this  study  do  not 
require  particular  values  in  any  axis,  the  two  groups 
of  data  were  considered  as  valid  separate  experiments 
and  evaluated  accordingly. 
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V.  Analysis  of  the  Data 

Statistical  Methods  Used  i a  Analysis 

The  first  step  in  analyzing  the  four  methods  was  to 
calculate  the  predicted  pilot  ratings  for  all  configur¬ 
ations.  These  were  then  compared  with  the  actual  pilot 
ratings  and  a  scatter  diagram  was  plotted  for  each 
method.  A  scatter  diagram  is  a  graphical  representation 
of  two  related  series  with  one  variable  scaled  on  the  X 
axis  and  the  other  scaled  on  the  Y  axis.  If  the  paired 
values  of  the  two  series  are  plotted  in  the  XY  plane, 
the  result  is  a  scatter  diagram.  In  this  thesis  the 
actual  pilot  ratings  were  scaled  on  the  X  axis  and  the 
predicted  pilot  ratings  were  scaled  on  the  Y  axis. 

A  line  of  perfect  correlation  was  then  plotted  on 
the  scatter  diagram.  This  line  would  be  the  locus  of 
all  points  if  there  were  perfect  correlation  or  agree¬ 
ment  between  the  predicted  and  actual  pilot  ratings. 
Since  both  axes  have  the  same  scale,  the  perfect 
correlation  curve  is  a  straight  line  represented  by  an 
equation  of  the  form 


where 


Y  =  a  +  bX 
a  =  0  and  b  =  1 . 
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Next,  a  line  of  regression  was  calculated  for  the 
plotted  data  points.  The  line  of  regression  is  a  curve, 
passing  through  the  plotted  points,  that  best  describes 
the  trend  of  the  data.  If  a  straight  line  trend  is 
assumed,  the  line  of  trend  will  be  represented  by  an 
equation  of  the  same  type  as  the  line  of  perfect 
correlation.  If  the  values  of  a  and  b  are  evaluated,  they 
can  be  compared  with  the"perfect"  values  of  the  correlation 
line.  To  solve  for  these  values,  the  principle  of  least 
squares  is  used.  This  principle  states  that  "A  line  of 
best  fit  to  a  series  of  values  is  a  line  the  sum  of  the 
squares  of  the  deviations  (the  differences  between  the 
line  and  the  actual  values)  about  which  will  be  a  min¬ 
imum"  (Ref  1:52). 

To  solve  for  a  and  b  the  following  values  are 
determined: 

«X(xy)  >  Y-W  '  N 

where  N  is  the  number  of  paired  values  being  considered. 
These  values  are  substituted  into  the  following  equations 
which  are  solved  simultaneously: 

YjX)  =  N  o  +  bYjx) 

^(XY)  =  a^TlX)  +  b^(X2) 
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The  line  of  regression  was  then  plotted  on  the 
scatter  diagrams.  Comparison  of  the  correlation  line 
and  the  regression  line,  which  will  be  discussed  more 
fully  in  subsequent  sections,  indicated  that  only  two 
of  the  methods  gave  good  results.  It  was  decided  to 
investigate  these  two  cases  further  by  calculating 
the  sample  variance  of  the  data  points  about  the  line 
of  perfect  correlation.  This  gives  a  numerical  value 
to  the  scatter  or  variation  of  the  data  points  about 
the  line  of  correlation,  which  is  a  figure  of  the 
relative  merit  of  the  methods  in  predicting  the  correct 
pilot  rating.  The  sample  variance  of  the  data  points 
about  the  line  of  perfect  correlation  was  calculated 
from  the  formula 


C Tc. 


■/«4  i 


a  4  I  o  r>) 


Analysis  of  Pilot  Consistency 

Figures  8,  9,  10,  and  11  show  the  pilot's  ability 
to  repeat  ratings  for  the  same  configurations.  The 
points  were  arrived  at  by  comparing  the  first  rating 
of  the  configuration  with  the  second,  the  second  with 
the  third,  and  so  on.  The  "old"  rating  was  scaled  on  the 
X  axis  and  the  "new"  rating  on  the  Y  axis.  Any  tendency 
of  the  pilot  to  change  his  rating  as  he  progressed  through 
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Single  Axis  Repetition  for  Three  Axis  Experiment 


Three  Axis  Repetition 
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the  sessions  is  reflected  in  the  shift  of  the  data 
points  above  or  below  the  line  of  perfect  correlation. 

An  examination  of  all  four  figures  shows  that  the 
pilot  had  a  tendency  to  rate  the  "poor"  configurations 
slightly  better  and  the  "good"  configurations  slightly 
worse  as  he  progressed  through  the  experiments.  The 
trend  of  the  ratings  is  seen  from  the  lines  of  regress¬ 
ion  which  are  plotted  as  solid  lines  on  all  the  figures. 
The  pilot  was  excellent  in  his  ability  to  repeat  ratings 
within  +1  rating  points.  All  of  the  cases  for  the  two 
axis  experiment  were  within  this  range.  In  the  three 
axis  experiment,  89%  of  the  single  axis  repeats  and 
91.5%  of  the  three  axis  repeats  were  within  +1. 

Analysis  of  Method  #1 

As  anticipated,  mtthod  #1  is  useless  for  predicting 
pilot  rating.  Figures  12,  13,  and  14  show  that  the 
line  of  regression  is  nowhere  near  the  perfect  correla¬ 
tion  line  and  the  method  consistently  predicts  ratings 
that  are  higher  than  the  actual  ones.  The  values  of 
(a)  and  <b)  for  the  line  of  regression  are  given  in 
Table  IV. 


Table  IV 

Method  #1,  Slope  (b)  and  Y  Intercept  (a) 
for  Line  of  Regression 


a 

b 

Two  axis  data 

4.84 

.6788 

Three  axis  data 

9.71 

.0325 

Combined  data 

6.32 

.4250 
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Analysis  of  Method  #2 

Method  92  is  one  of  the  two  methods  that  is  con¬ 
sidered  promising.  Although  the  regression  line  for 
the  two  axis  experiment  (Pigure  IS)  is  skewed  from 
the  correlation  line,  the  data  is  well  grouped  around 
it.  In  the  three  axis  cases  (Pigure  16)  there  is  a 
definite  trend  to  predict  high  by  .8  of  a  rating  point. 
However,  if  the  two  and  three  axis  data  is  combined 
(Pigure  17)  to  increase  the  range  of  points,  the  trend 
and  correlation  line  closely  approximate  each  other. 
There  is  a  gradual  tendency  to  predict  high  as  the  con¬ 
figurations  become  poorer.  The  values  for  (a)  and  (b) 
for  method  92  are  given  in  Table  V. 


Table  V 

Method  92,  Slope  (b)  and  Y  Intercept  (a) 
for  Line  of  Regression 

a 

b 

Two  axis  data 

2.00 

.5997 

Three  axis  data 

.77 

1.0011 

Combined  data 

.03 

1.0832 

Analysis  of  Method  9 3 

Method  93  is  the  second  method  that  appears  to  be 
promising.  The  predicted  ratings  for  the  scatter 
diagrams  were  determined  from  the  curve  of  Pigure  C-l, 
which  is  the  approximate  center-line  of  Ashkenas'  band 
curve  as  taken  from  reference  11. 


34 


GS/EB/62-4 

The  trend  of  the  points  is  reasonably  close  to 
the  correlation  line  over  the  range  of  the  data  for 
the  two  and  three  axis  experiments  (Figures  18  and  19). 
Combining  the  two  and  three  axis  data  to  extend  the 
range  gives  a  line  of  regression  that  agrees  almost 
perfectly  with  the  correlation  line  (Figure  20).  The 
values  for  (a)  and  (b)  for  method  #3  are  given  in  Table 
VI. 


Table  VI 

Method  #3,  Slope  (b)  and  Y  Intercept  (a) 
for  Line  of  Regression 


a 

b 

Two  axis  data 

1.54 

.8142 

Three  axis  data 

-4.05 

1.5028 

Combined  data 

.17 

1.0167 

Analysis  of  Method  £4 


The  scatter  diagrams  indicate  that  method  #4  tends 
to  predict  higher  than  actual  ratings.  Ihe  trend  is 
not  as  pronounced  as  in  method  #1,  but  even  combining 
the  data  (Figure  23)  from  the  two  and  three  axis  cases 
(Figures  21  and  22)  produces  a  regression  line  that 
indicates  a  consistent  prediction  of  one  rating  point 
high.  This  bias  might  be  attributed  to  a  poor  choice 
of  the  "best"  rating.  However,  in  this  study  extra  runs 
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in  this  configuration  were  flown  to  insure  the  "best" 
cases  were  well  established.  The  values  for  (a)  and 
(b)  are  given  in  Table  VII. 


Table 

Method  #4,  Slope  (b) 
for  Line  of 

VII 

and  Y  Intercept 
Regression 

(a) 

a 

b 

Two  axis  data 

1.67 

.7887 

Three  axis  data 

4.14 

.6456 

Combined  data 

1.19 

.9788 

Comparison  of  Methods  #2  and  #3 

Because  of  the  paucity  of  data  in  the  two  axis 
experiment  and  the  limited  range  of  data  in  the  three 
axis  experiment,  there  is  not  sufficient  justification 
at  this  time  to  make  a  comparison  as  to  which  method 
would  be  better  for  either  case.  However,  the  combined 
data  is  compared  with  the  perfect  correlation  line  for 
values  of  (a),  (b),  and  sample  variance  in  Table  VIII. 

The  comparison  shows  that  both  methods  are  very  good 
for  predicting  the  correct  rating  and  there  is  very 
little  difference  between  the  scatter  of  their  points 
about  the  line  of  correlation. 

Compared  on  the  basis  of  ease  and  speed  of  calculation, 
method  #2  is  better.  However,  it  has  a  major  disadvantage, 
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since  a  "best”  combined  axis  pilot  rating  must  be 
determined.  On  the  other  hand,  method  #3  requires 
more  time  for  calculations,  but  can  be  used  if  only 
single  axis  data  is  available.  Because  most  of  the 
experimental  data  in  the  handling  qualities  field  has 
been  for  single  axis  cases,  method  #3  is  probably  the 
most  useful. 


Table  VIII 

Comparison  of  Line  of  Perfect  Correlation 
Line  of  Regression  and  Sample  Variance 
of  Method  #2  and  Method  #3 

with 

a 

b 

Line  of  perfect  correlation 

0 

1.0 

0 

Method  it  2 

.03 

1.0832 

.8838 

Method  ft 3 

.17 

1.0167 

.8469 
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VI .  Conclusion  and  Recommendations 

Conclusion 

This  thesis  has. proven  conclusively  that  a  predicted 
pilot  rating  for  a  multiple  axis  system  cannot  be  found 
by  adding  the  pilot  ratings  of  the  individual  axes  used 
in  the  system  (method  #1). 

It  has  shown  that  there  are  at  least  two  methods  for 
determining  a  predicted  rating  that  give  good  results. 
These  two  methods  (methods  #2  and  #3)  were  compared  and 
it  was  concluded  that  method  #3  was  probably  the  most 
useful,  due  to  the  type  of  data  that  could  b-^  used  with 
it. 

It  has  shown  that  a  fourth  method,  which  utilized 
features  of  methods  #2  and  tt 3,  gave  poorer  results  than 
the  individual  methods. 

Recommendations 

A  follow-on  experiment  is  recommended  to  collect  data 
on  more  test  subjects.  Investigation  could  be  concen¬ 
trated  on  the  two  most  promising  methods  to  determine  if 
one  gives  consistently  better  results  for  a  larger  sample 
of  pilots.  Techniques  to  determine  the  "best"  combined 
axis  rating  used  in  method  #2  from  single  axis  data  could 
be  investigated.  Additional  studies  to  determine  if  one 
method  gives  better  results  as  the  number  of  axes  are 
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changed  could  be  made.  The  conversion  curve  used  in 
method  #3  could  be  refined  to  determine  if  the  center- 
line  plot  used  in  this  thesis  gives  the  best  results. 
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Appendix  A 

Desc  rintion  of  Expo  rlmental  Equipment 

The  equipment  used  during  the  experiment  consisted 
of  two  analog  computers,  a  cockpit  mock-up  in  which  a 
sidestick  and  rudder  pedals  were  installed,  an 
oscilloscope,  a  galvanic  type  meter  (sideslip  indicator), 
two  low  frequency  function  generators,  an  audio 
oscillator,  a  six  channel  oscillograph,  and  simple 
passive  circuit  elements.  A  general  view  of  the  experi¬ 
mental  set-up  is  shown  in  Figure  A-l. 

Analog  Computer 

The  control  and  reference  voltage  circuits  of  the 
two  Goodyear  Model  L-3  analog  computers  were  connected 
to  operate  both  units  from  a  single  control.  Input 
generators  #3  and  #4  in  Figure  A-2  were  wired  on  one 
computer  anc  provided  sine  and  cosine  signals  with 
amplitudes  and  frequencies  of  1.15  volts  at  .816 
radians/second  and  1.64  volts  at  1.32  radians/second. 

The  controlled  elements  J^- ,  £  »  the  sum¬ 

ming  points  were  mechanized  using  the  rest  of  the  com¬ 
puter  elements.  Schematic  diagrams  of  and 

1  are  shown  in  Figures  A-3,  A-4,  and  A-5.  The  values 
<5"> 

for  -  — ,  f  ,  and  f  were  varied  with  pots  A,  B,  C,  and 
D  of  the  schematics.  Figure  A-2  shows  the  summing  points 
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used  to  generate  the  random  input  signals.  Input 
generators  #1  and  #2  were  Hewlett-Packard  function 
generators  producing  sine  waves  of  4.8  volts  at  .367 
radians/second  and  5.5  volts  at  .483  radians/second. 

The  inputs  were  compared  with  the  outputs  of  the  con¬ 
trolled  elements  (©,<£,  $  )  and  the  results  (  ©c , 

<Pf  t  &  )  were  fed  to  the  displays. 

Because  of  the  sensitivity  of  the  sideslip  needle, 
the  £  voltage  was  reduced  by  a  40kilohm  resistor  which 
brought  the  needle  deflections  within  the  desired  range. 
The  pitch  error  voltage  was  reduced  through  a  sign 
changer  and  summed  with  the  output  of  the  bank  error 
circuit.  The  operation  of  the  bank  circuit  is  discussed 
in  detail  in  the  next  section. 

Pitch  and  Bank  Display  Circuit 


A  realistic  pitch  and  bank  display  was  generated  by 
the  circuit  shown  in  Pigure  A-6. 
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A  50  cycle  sine  wave  was  fed  to  the  horizontal 
sweep  of  the  oscilloscope  and  to  an  electronic  multi¬ 
plier  where  it  was  multiplied  by  the  bank  error  signal, 

<P£  .  The  results  of  the  multiplication  was  then 
summed  with  the  pitch  error  signal,  ©f  ,  and  fed  to 
the  vertical  input  of  the  scope.  Because  of  the  phase 
relationships  of  the  horizontal  and  vertical  inputs, 
the  bank  angle  of  the  bar  was  a  direct  function  of  the 
bank  error  signal.  The  vertical  movement  of  the  center 
of  the  bar  was  a  direct  function  of  the  pitch  error 
signal.  The  HC  circuit  was  placed  in  the  input  lead 
of  the  50  cycle  signal  to  the  electronic  multiplier  to 
compensate  for  the  phase  shift  introduced  by  the  multi¬ 
plier. 

Cockpit  Mock-up 

The  location  of  the  displays  and  control  devices  in 
the  cockpit  are  shown  in  Figures  A-7  and  A-8.  A  detailed 
sketch  of  the  sidestick  is  illustrated  in  Figure  A-9. 

Control  Devices 

The  pilot  controlled  pitch  and  bank  with  a  side-stick 
mounted  as  illustrated  in  Figure  A-7.  The  sidestick,  an 
B-6  Autopilot  Control  Stick,  manufactured  by  Minneapolis- 
Honeywell,  had  been  modified  by  removing  rate  damping 
devices,  leaving  only  spring  restrainers.  The  maximum 
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deflection  of  +20°  in  each  axis  required  approximately 
2%  pounds  pressure  at  the  handgrip.  Because  of  the 
greater  distance  to  the  axis  of  lateral  movement  than 
to  the  axis  of  longitudinal  movement  the  pilot  was 
required  to  move  the  stick  a  greater  linear  distance 
laterally  than  longitudinally  for  the  same  amount  of 
angular  deflection.  The  stick  had  a  mechanical  dead¬ 
band  of  approximately  +  1°  longitudinally  and  + 
laterally. 

The  rudder  pedals,  used  to  control  the  sideslip 
needle,  were  spring  restrained  as  shown  in  Figure  A-10. 
A  deflection  of  approximately  +_  2  inches  required  about 
40  pounds  of  pressure.  The  rudder  force  versus  rudder 
deflection  characteristics  are  shown  in  Figure  A-ll. 

The  wiring  for  the  control  devices  is  shown  in 
Figure  A-12.  Bach  axis  of  the  sidestick  was  mechanic¬ 
ally  linked  to  the  arm  of  a  400  ohm  potentiometer.  The 
rudder  pedals  were  connected  mechanically  through  a 
pulley  arrangement  to  the  arm  of  a  50  kilohm  rotary 
potentiometer.  The  electrical  outputs  of  the  potentio¬ 
meter  arms  were  fed  to  amplifiers  1,  2,  and  3  (Figures 
A- 3,  A- 4,  and  A-5).  Three  trim  potentiometers,  tied 
in  parallel  with  the  control  potentiometers,  were  also 
fed  to  these  amplifiers  to  provide  a  means  of  bringing 
the  electrical  input  to  zero  when  the  control  was  not 
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deflected.  Power  to  the  potentiometers  was  supplied 
by  the  computer  reference  voltages  (+200  VDC).  How¬ 
ever,  because  of  the  low  power  ratings  of  the  400  ohm 
pots,  it  was  necessary  to  reduce  this  to  +5  volts  with 
large  variable  resistors  in  series.  The  computer 
ground  was  used  as  center  tap  for  the  pots  with  the 
opposite  polarities  applied  to  each  end.  Two  Zener 
diodes  kept  the  voltage  regulated  very  close  to  5 
volts.  The  voltage  output  versus  physical  deflection 
characteristics  of  the  control  devices  illustrated  in 
Figures  A-13,  A-14,  and  A-15,  were  measured  at  the  out¬ 
put  leads  of  Figure  A-12  with  a  constant  +5  VDC  across 
the  pots. 


System  Sensitivity 

The  display  sensitivities  were 
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The  overall  sensitivities  from  control  device  to 
display  were 
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General  View  of  Simulation  Equipment 
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Schematic  Diagram  of  Simulated  Controlled  Element 
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Schematic  Diagram  of  Simulated  Controlled  Element 
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Figure  A- 5 

Schematic  Diagram  of  Simulated  Controlled  Element 
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Close-up  View  of  Cockpit 


Figure  A-8 

Sketch  of  Cockpit  Equipment  Layout 
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Close-up  View  of  Rudder  Pedals 
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Rudder  Deflection  versus  Applied  Pressure 
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Schematic  Diagram  of  Wiring  for  Control  Devices 
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Sidestick  Deflection  versus  Potentiometer  Output 
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Left-Right  Sidestick  Deflection  versus  Potentiometer  Output 
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Rudder  Deflection  versus  Potentiometer  Output 
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Appendix  B 

Determination  of  System  Dynamics 

Derivation  of  Transfer  Functions 

The  transfer  functions  used  in  this  experiment  were 
derived  from  the  following  equations  of  motion  for  a 
tandem  rotor  helicopter: 

Longitudinal 

xi  +  32.2  O  =  -\0  4  SxX.  -+-.0 1 1/x)  1. 4  8  Q  4-  .463  8  es 

x^'2310  =  -.300^ — 1.816+. 494 

a  ='-.OC385yav.OH8^r-l.65©^.53i  <5 

5 

Lateral  -  Directional 

» 

^'32.2^231  ^  =  -1  52/tr-2.  55<p-b/.Z  8a, +  .0*2 Sr? 

£-.671  ^  =  - 00855  /zr-i  260V-.OI87  ^647^  “2  2 1  <5r, 
^-.088£=y-.OO374/3r  -.0084^-035  .0057 S^+.ISZ  6<r 
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Substituting  /trw U0q  and  solving  for  the  transfer  functions 
gave  the  following  results: 


6S)  -  K(S+.Q4  2  6)  C 5  4-.9U) _ 

&.5fe)  (S  -.628)  (S  + 127)  [?t2O6^20iThi09*I 

0CS)  _  K(S  -630(Sf  0 _ 

£,  (S)  "  (S+.549KS  +  l.?84)(S-.OI6)fS-7M) 


$(S)  _  K  (S+,00 1  2  9~)  (S  -H.4  |)(S  -  1 3S5) 

A  ts>  -  (S-.0 1 6)  (S-.784)(  S+.S  49)  (S  +  1.7  84) 


To  determine  what  effect  augmentation  had  on  the 
poles  and  zeros  of  the  transfer  functions,  rate  feed¬ 
back  in  pitch  and  yaw  were  introduced  into  the  control 
loop  (Figure  B-l). 

The  sensitivity,  SA—  is  the  sensitivity  of  the 
controlled  part  of  the  airframe  to  the  movement  of  the 
control  device  and  is  considered  to  be  a  pure  gain  term. 
Then  in  general 


&  =  SA[Sj  -  Ka#s  Ao]  _  SA[8j  -  K^Ao] 
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Closed-loop  Syste:.:  Using  Rate  Feedback 
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and 


&e  S©  [  S  es  "  K®  0  ] 

6r  =  S*  [  brf  ~ 

Substituting  Se  and  S*.  into  the  general  form  of  the 
equations  of  motion  for  this  aircraft 


sOi  -4  cj  ©  '  ^  AX  M-  4  X  SAJ-  A/J~  +  iff  &  *  ^$e  “  K*  ©3 

Alt'  ~  Oo  ©  ~  ^aa  AJ-  ~*~Z  aaT  sU)  4-  &  +Z 

©  +M*  Q  Ki©)l 

>zk- +■  Uo^ 31  vi +y+$  +  t  wh tv csrp~K^  t)] 

<f>  -Z_£*  f  -  L/v-AT'-t  L i (j>  +  V  t L^U+ L^tV'fir,-  J 

lx* 

P-hi  4  -  A4r  ^  *-/V^'  ^*A/f  ^-#-A/fc  &  +  M*rLV&r,-fcf  ft] 
r«? 
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Rearranging  terms 


yCA.  V  0  =  Xu  M-  +  +  (Xi-Se  X^t  i^e)©  +•  Xje  <§e, 

-  Uo  &  ~£u.  AX  +2mtM~+  -  Se  2^,  Ke)©t  Se;?s«  Se, 

e  t  (Hi  -  S. i  Hi.  Ki)  ©  *  5 ,M{. 

/i) - <^f>  -(a,  *r  -tVpjt  *  ^(.-.(  U.tSrV(r.tCi)'t'-f-Sr  HrSff 

$  ~  llLS.  Lsr  ArfLfp  j>  +  L £  Sa  +  (L+  -  <5V^jr  &)  ^  <5V  ^$r  <SV^ 

lx* 

-  £*Z.  cf)  -^4  ^fhJjf  j>  +  A/fr  Sa+(N+  ~  Sfhllr  +  S+Nzr’&rf 

£  it 

Referring  to  the  actual  equations  of  motion 


Se  X,.=  .463 

and 

*•'  = 

“0.4  8 +  .4  6  3  K*) 

Se  H„  =  .494 

and 

z;  = 

“0.5  1  +4  94  K*) 

—  •'531 

and 

m;= 

“0.85  +  .5  31  KJ 

SY  Yw  =  .042 

and 

u;  = 

(2  31  +  .04  2  K,) 

S*  =-.22  1 

and 

/ 

L*.  = 

CO  1  87+2  2  1  KJ 

S„  Ntr=  .2  52 

and 

n'  = 

—  (035  +  .2  5  2  K*) 
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To  establish  an  initial  value  for  K*  the  short 
period  frequency  was  arbitrarily  set  equal  to  zero. 

Using  the  short  period  approximation  (Ref  8:111-7) 

u>.Pa  =  <  -Ue  M_=0 

and 

M*'  =  Ug-M^:.  =  -3.66  5 

Substituting  into  the  equation  for  M 

K6i=  3.4  2  5 

Multiples  of  the  initial  K*  were  then  used  to  find 
additional  values  of  mJ  which  were  substituted  into 
the  original  equations  of  motion.  The  transfer  functions 
determined  are  listed  in  Table  B-I. 

Analysis  shows  the  positive  pole  moving  toward  the 
origin  as  Kj  is  increased.  This  pole  was  used  to  vary 
the  amount  of  instability  in  pitch.  Since  it  was  not 
necessary  to  simulate  a  particular  aircraft,  the  other 
poles  and  zeros  were  arbitrarily  fixed  at  values 
approximating  those  for  K©  =  3.425.  The  transfer 
function  used  in  pitch  was 
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9(5)  _  K.rs  +  .04)(S+.9) _ 

6-.CS)  _  (S  +  U  )  (S  +  5)  [Si+2C7)C2S)S  +.2  S'J 

•p. 


Table  B-I 

Transfer  Function  f for  Different  K* 

Qtf(S) 


“K*  Transfer  Function 


1V  • 

0 

K(S+. 0426 )(S+. 914) 

(S- . 629CS+3.27 )  [Sa  +2(  .  364) (  . 209 )S+ . 2094] 

K( S+. 0426 )(S+. 914) 

3.425 

(S-.287)(S+4.58)  [S*+2(  .611)(.261)S  +  .261*J 

13.7 

K(S+. 0426 )(S+. 914) 

(S-.lll)(S+9. 51)(S+. 504)(S+.168) 

27  4 

K( S+. 0426 )(S+. 914) 

(S-.0689)(S+16.6)(S+.696)(S+.113) 

A  similar  approach  was  used  in  determining  values 
for  Lv'  and  N*.'  ,  but  K*  was  arbitrarily  given  values 
of  5,  10,  and  20.  The  transfer  functions  determined 
are  listed  in  Tables  B-n  and  B-III. 
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Table  B-II 


Transfer  Function  for  Different  Ky 

K*  Transfer  Function 


K(S-.836)(S+1) 

0 

(S-.016)(S-.784)(S+.549)(S+1.784) 

K(S-.37)(S+1.8) 

5 

(S+1.15XS+2.18)  LSl  ~2(  . 546 )(  .  462)S+  .462'1] 

K(S-.179)(S+2.87) 

10 

(S+1.38XS+2.79)  [.S1  - (  .  340 ) (  .  495 )S+ .  495*] 

K(S-.0359)(S+5.25) 

20 

(S+1.51XS+5.18)  [S1-2(.205)(.513)S  +  .513*] 

As  K*  increases  the  set  of  complex  positive  poles 
moves  to  the  left.  The  stability  of  this  axis  was 
varied  by  changing  the  damping  ratio  of  these  poles. 

A  further  simplification  was  made  by  cancelling  the 
pole  and  zero  that  are  approximately  equal  for 
Ky  =10  and  20.  The  positive  zero  was  moved  to  the  left- 
haid  side  of  the  s-plane  to  improve  the  overall  closed 
loop  stability.  The  final  form  of  the  bank  equation  was 


<£(S)  _  Ka(S  -U) _ 

d,(S)  (s  +  1.5)  [SX+2%C5)S  +.5*} 
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Table  B-III 

Transfer  Function  f for  Different  K* 

o  r-(0 


Transfer  Function 


0 

K(S+. 00129 )(S+1. 41) (S-1355) 

(S-.016)(S-.784)(S+.549)(S+1.784) 

K(S+.  00129  )(S+1. 41)  (S-1355) 

5 

<S+1.15)(S+2.12)  LS‘  -2(.546)(.462)S+.462I] 

K.(S+  .0220)  (S+l.  30)  (S- 1353) 

10 

(S+1.30)(S+2.97)  ]S‘-2(.340)(.495)S+.  495*1 

K(S+.001)(S+1.41)(S-1351) 

20 

(S+1.51)(S+5.18)  [S  ^ -2(.  205  )(.  513  )S+.  513*1 

The  same  procedure  was  used  in  determing  that  the  damp¬ 
ing  ratio  would  be  used  in  adjusting  the  stability  in 
sideslip  and  the  pole-zero  set  was  again  cancelled.  The 
large  positive  zero  has  the  effect  of  making  the  function 
negative.  The  function 

PCS)  _  K.  5 _ 

6>(s)  (s  +  Sfs'+^csis  +.sJ] 


was  adopted  for  the  initial  trials  in  the  simulator. 
However,  the  operator  was  not  able  to  control  the  side¬ 
slip  indicator  within  desirable  display  limits.  There¬ 
fore,  the  function 

$(5>_  K.S  -*•  •>) _ 

&rp(S)  (S  +  3)[S'  +  ^-f,(5)S  +  .5’] 
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Determination  of  Augmentation  Levels 

The  particular  settings  on  the  computer  for  each 
augmentation  level  were  selected  by  first  estimating 
values  of  ~ —  ,  f.  ,  and  /  that  would  produce  a 

TC# 

rating  in  each  major  category  of  the  Cooper  scale.  A 
pilot  (the  author)  who  was  not  flying  the  data  runs 
then  flew  and  rated  single  axis  runs  for  these  settings. 
Adjustments  were  made  in  the  settings  until  he  was 
rating  the  "H"  level  at  2,  and  "M"  level  at  4t  and  the 
"L"  level  at  6.  The  values  at  these  settings  were  then 
fixed  for  the  rest  of  the  experiment. 
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Appendix  C 

Data  and  Curve  Used  in  Calculations 
Conve rsion  Curve 

The  curve  used  to  determine  the  predicted  pilot 
ratings  in  method  #3  and  #4  is  shown  in  Figure  C-l. 

It  is  the  approximate  center  line  of  the  band  curve 
proposed  by  Ashkenas  in  reference  11. 

Raw  Data 

The  raw  data  taken  during  the  experiment  is  listed 
in  Table  C-I.  It  is  available  for  checking  the  cal¬ 
culated  results  shown  in  this  thesis.  It  can  also  be 
used  for  other  problems. 

Numbers  Used  In  Calculations 

The  numbers  used  in  all  calculations  were  averaged 
values  of  the  ratings  for  any  particular  configuration. 
These  numbers  are  given  in  Tables  C-II,  C-III,  C-IV,  and 
C-V. 
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ftin  Configuration  Pilot  Rating 


©  <P  £ 


Data  Run  #2 

Three  Axis  -  120  seconds 
ftrn  Configuration  Pilot  Ratini 

e  <P  g _ 

22  H  L  H  -* 

23  H  L  H  9.0 
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Data  Runs 


Run  Configuration  Pilot  Ratinj 


& 

<*> 

& 

12 

H 

3.0 

13 

L 

7.0 

14 

H 

2.5 

15 

M 

4.5 

16 

M 

3.5 

17 

L 

5.5 

18 

H 

3.0 

19 

H 

3.5 

20 

L 

7.0 

21 

M 

3.0 

Rin  Configuration  Pilot  Ratini 

Q  <fi  g _ 

24  M  M  M 

25  M  M  M  8.0 


Pilot  exercised  re-run  option 


GE/E E/62-4 

Run  Configuration  Pilot  Rating  foin  Configuration  Pilot  Rating 


© 

<P 

f 

© 

P 

e 

26 

L 

M 

M 

- 

49 

L 

H 

M 

10.0 

27 

L 

M 

M 

9.0 

50 

M 

H 

H 

6.5 

28 

H 

H 

L 

- 

51 

H 

H 

M 

6.5 

29 

H 

H 

L 

10.0 

52 

L 

L 

L 

- 

30 

L 

L 

M 

- 

53 

L 

L 

L 

- 

31 

L 

L 

M 

10.0 

54 

L 

L 

L 

10.0 

32 

H 

L 

M 

- 

55 

H 

M 

H 

6.0 

33 

H 

L 

M 

8.0 

56 

L 

L 

H 

- 

34 

M 

L 

L 

- 

57 

L 

L 

H 

9.5 

35 

M 

L 

L 

10.0 

58 

L 

H 

H 

- 

36 

H 

M 

L 

- 

59 

L 

H 

H 

10.0 

37 

H 

M 

L 

9.0 

60 

H 

L 

L 

- 

38 

H 

H 

H 

6.0 

61 

H 

L 

L 

9.5 

39 

L 

M 

H 

- 

62 

M 

L 

M 

«■ 

40 

L 

M 

H 

9.0 

63 

M 

L 

M 

9.5 

41 

M 

H 

M 

8.0 

64 

L 

H 

L 

- 

42 

M 

M 

H 

7.5 

65 

L 

H 

L 

10.0 

43 

L 

M 

L 

- 

66 

M 

M 

L 

mm 

44 

L 

M 

L 

9.5 

67 

M 

M 

L 

9.0 

45 

M 

H 

L 

- 

68 

H 

M 

M 

7.0 

46 

M 

H 

L 

9.5 

69 

M 

L 

H 

mm 

47 

L 

H 

M 

- 

70 

M 

L 

H 

8.5 

48 

L 

H 

M 

- 
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H 

H 

H 

5.5 
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Data  ftin  #3 

Three  Axis  Repeated  Runs  -  120  seconds 

Run  Configuration  Pilot  Rating  Run  Configuration  Pilot  Rating 


e 

$ 

e 

4> 

e 

72 

H 

L 

H 

- 

90 

H 

H 

L 

9.0 

73 

H 

L 

H 

7.0 

91 

L 

L 

H 

- 

74 

L 

H 

H 

- 

92 

L 

L 

H 

8.5 

75 

L 

H 

H 

9.0 

93 

M 

L 

H 

8.0 

76 

L 

M 

L 

- 

94 

H 

M 

M 

6.0 

77 

L 

M 

L 

9.0 

95 

H 

H 

L 

- 

78 

M 

H 

H 

6.5 

96 

H 

H 

L 

10.0 

79 

L 

L 

M 

- 

97 

M 

H 

L 

- 

80 

L 

L 

M 

- 

98 

M 

H 

L 

- 

81 

L 

L 

M 

9.0 

99 

M 

H 

L 

9.0 

82 

M 

H 

M 

7.5 

100 

H 

M 

L 

- 

83 

H 

H 

H 

6.0 

101 

H 

M 

L 

w* 

84 

M 

M 

M 

7.0 

102 

H 

M 

L 

8.0 

85 

H 

H 

M 

6.0 

103 

H 

M 

H 

5.0 

86 

L 

H 

M 

- 

104 

M 

H 

M 

6.5 

87 

L 

H 

M 

8.5 

105 

M 

H 

H 

5.5 

88 

H 

H 

L 

- 

106 

H 

M 

H 

5.5 

89 

H 

H 

L 

— 

107 

H 

H 

H 

6.0 
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Data  Ain  #4 

Single  Axis  Repeated  Runs  -  60  seconds 

Run  Configuration  Pilot  Rating  Ain  Configuration  Pilot  Rating 


e 

<t>  $ 

© 

e 

108 

M 

4.0 

120 

L 

6.0 

109 

L 

5.0 

121 

H 

3.0 

110 

M 

4.5 

122 

M 

5.0 

111 

H 

3.5 

123 

H 

3.0 

112 

M 

4.5 

124 

H 

4.0 

113 

L 

6.0 

125 

M 

4.5 

114 

L 

6.0 

126 

M 

5.0 

115 

M 

4.5 

127 

L 

6.5 

116 

M 

5.0 

128 

L 

5.0 

117 

H 

4.0 

129 

H 

3.5 

118 

H 

3.0 

130 

H 

3.0 

119 

M 

5.0 

131 

M 

4.0 

Two 

Axis 

Experiment 

Data  Run  #5 

Single  Axis  -  60  seconds 

Ain 

Configuration  Pilot  Rating 

Ain 

Configuration  Pilot  Rating 

© 

<t> 

© 

132 

H 

3.0 

136 

H 

3.0 

133 

M 

3.5 

137 

M 

3.5 

134 

L 

4.0 

138 

H 

3.0 

135 

L 

3.5 

139 

H 

2.5 

80 


GE/EE/62-4 

Run  Configuration  Pilot  Rating  Ran  Configuration  Pilot  Rating 


© 

4> 

© 

140 

L 

4.0 

145 

H 

3.0 

141 

H 

3.0 

146 

H 

4.0 

142 

M 

4.5 

147 

L 

6.0 

143 

M 

5*0 

148 

L 

5.5 

144 

L 

6.5 

149 

L 

5.0 

Data  Rin 

#6 

2  Axis  - 

120 

seconds 

Run 

Configuration  Pilot  Rating 

Rin 

Configuration 

Pilot  Rating 

e 

<P 

© 

4> 

150 

M 

M 

4.0 

165 

L 

L 

7.0 

151 

H 

M 

3.0 

166 

H 

M 

4.0 

152 

H 

H 

2.0 

167 

H 

H 

3.0 

153 

M 

H 

3.5 

168 

M 

M 

4.0 

154 

L 

H 

4.0 

169 

M 

M 

5.0 

155 

M 

L 

6.5 

170 

L 

L 

7.0 

156 

L 

M 

4.0 

171 

M 

H 

3.0 

157 

L 

L 

7.0 

172 

H 

H 

3.5 

158 

H 

L 

6.0 

173 

L 

M 

4.0 

159 

H 

M 

4.0 

174 

L 

M 

5.0 

160 

H 

H 

3.0 

175 

H 

L 

6.5 

161 

M 

H 

3.0 

176 

L 

H 

4.0 

162 

L 

H 

4.0 

177 

M 

L 

7.0 

163 

H 

H 

3.0 

178 

H 

L 

6.0 

164 

M 

L 

7.0 
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Table  C-II 

Averaged  Values  of  Single  Axis  Runs 
for  Two  Axis  Experiment 

Pitch  Bank 


H 

-  2.9 

H  -  3.3 

M 

-  3.5 

M  -  4.8 

L 

-  3.8 

L  -  5.8 

Table  C-III 

Averaged  Values  for  Two  Axis  Runs 

Configuration 

Pilot  Rating  Configuration 

Pilot  Rating 

© 

© 

H 

H 

2.9  M  L 

6.8 

H 

M 

3.7  L  H 

4.0 

H 

L 

6.2  L  M 

4.3 

M 

H 

3.2  L  L 

7.0 

M 

M 

4.3 

Table  C-IV 

Averaged  Values  of  Single  Axis  Rms 
for  Three  Axis  Experiment 

Pitch 

Bank  Sideslip 

H 

-  2.8 

H  -  3.3 

H  -  3.4 

N 

-  4.6 

M  -  4.3 

M  -  4.1 

L 

-  6.7 

L  -  5.8 

L  -  6.2 
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Table  C-V 

Averaged.  Values  of  Three  Axis  tons 
Configuration  Pilot  Rating  Configuration 


e  <t>  $ 

H  H  H 
H  H  M 
H  H  L 
H  M  H 
H  M  M 
H  M  L 
H  L  H 
H  L  M 
H  L  L 
M  H  H 
M  H  M 
M  H  L 
M  M  H 
M  M  M 


5.9 

6.2 

9.7 

5.5 

6.5 

8.5 
8.0 
8.0 

9.5 
6.2 
7.3 
9.2 

7.5 
7.5 


©  ^  ^ 

M  M  L 

M  L  H 
M  L  M 
M  L  L 
L  H  H 
L  H  M 
L  H  L 
L  M  H 
L  M  M 
L  M  L 
L  L  H 
L  L  M 
L  L  L 


Pilot  Rating 


9.0 

8.2 

9.5 

10.0 

9.5 

9.2 

10.0 

9.0 

9.0 

9.2 

9.0 

9.5 

10.0 
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